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Abstract 


An investigation into the aerodynamics of the Launch Abort Vehicle for NASA’s Constel- 
lation Crew Launch Vehicle in the subsonic, incompressible flow regime was conducted in 
the NASA Langley 20-ft Vertical Spin Tunnel. Time histories of center of mass position 
and Euler Angles are captured using photogrammetry. Time histories of the wind tunnel’s 
airspeed and dynamic pressure are recorded as well. The primary objective of the inves- 
tigation is to determine models for the aerodynamic yaw and pitch moments that provide 
insight into the static and dynamic stability of the vehicle. System IDentification Programs 
for AirCraft is used to determine the aerodynamic model structure and estimate model pa- 
rameters. Aerodynamic models for the aerodynamic body Y and Z force coefficients, and 
the pitching and yawing moment coefficients were identified. 
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1 Introduction 


In 2004, President Bush [1] outlined a plan for the future exploration of outer space. 
One of the major objectives of this plan is to return to the moon by 2020. In response to 
this vision, NASA started to design and develop a new launch system to replace the space 
shuttle. The architecture for this system is defined in the Exploration Systems Architecture 
Study (ESAS) final report. [2] The ESAS final report recommended that the new launch 
system needs a Launch Abort System (LAS) capable of pulling the Crew Exploration Vehicle 
(CEV)’s Crew Module away from the Crew Launch Vehicle (CLV) during launch failures. 
While the Launch Abort Vehicle (LAV) (See Figure 1) being developed will incorporate 
many modern technologies, the concept is similar to the Launch Escape System (LES) 
utilized during the Mercury and Apollo Programs. [3, 4] 



Figure 1. General configuration of Launch Abort Vehicle. 

An investigation into the low subsonic aerodynamic stability characteristics of the NASA 
Orion Launch Abort Vehicle (LAV) was conducted in NASA Langley Research Center 
(LaRC) 20-foot VST. The goal of this investigation is to model the aerodynamic forces 
and moments acting on the LAV in the subsonic, incompressible flow regime. The primary 
objective is to determine models that describe the longitudinal and directional stability of 
the LAV. A secondary objective is to determine models for the remaining aerodynamic 
forces and moments. 

This investigation focuses on the low-speed (Mach 0.2) heat shield-forward flight of the 
LAV during a pad abort scenario. During this scenario, the LAV is thrust away from the 
CLV by an abort motor located in the LAS. The LAV will be flying tower-end forward 
during this time. The LAV must reorient to a stable heat shield-forward flight state so that 
the LAS can be jettisoned from the Crew Module to allow for a safe return to the ground. 
Understanding the LAV’s aerodynamic characteristics after reorientation and before LAS 
jettison when the LAV is at low altitude and speed is critical to ensuring a safely controlled 
descent. 

This investigation was carried out in the NASA LaRC 20-ft VST. A dynamically scaled 
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model of the LAV was allowed to fly unconstrained inside the wind tunnel. Data analysis 
was carried out using System IDentification Programs for AirCraft (SIDPAC) [5]. SIDPAC 
is a collection of software tools developed by Dr. Vladisav Klein and Dr. Eugene Morrelli 
at National Aeronautics and Space Administration (NASA) LaRC for conducting system 
identification. The collection consists of different output-error and equation-error estimation 
algorithms developed for both the time and frequency domains. Additional tools for deter- 
mining derivatives, Fourier transforms, and data smoothing are also included in SIDPAC. 
While some of the SIDPAC algorithms are written specifically for aircraft, many of the 
core algorithms were developed in a generalized form which allows flexability in their im- 
plementation. The analysis presented in this report demonstrates how SIDPAC’s analysis 
capabilities can be extended beyond aircraft. 

This report documents the approach used to determine the aerodynamic force and mo- 
ment models. Section 1 describes the scope of the investigation and outlines the report. 
Section 3 discusses the method used to analyze the collected data. Section 2 discusses the 
facilities and equipment used to obtain data. Section 4 provides the results of the data 
analysis and discusses the meaning of the results. Section 5 provides a concluding summary 
of the report. 
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2 Experiment 

2.1 Test Facilities 

The experimental portion of this investigation took place at the NASA LaRC’s 20- foot 
Vertical Spin Tunnel (VST) [6]. The VST has been in use by NASA LaRC since 1941. It 
has been used to analyze the spin and spin recovery characteristics of many aircraft. This 
facility has participated in investigations for NASA’s Mercury [7], Gemini [8], Apollo [9], 
Stardust [10-12], and Mars Sample Return [13] missions. Basic research concerning subsonic 
blunt-bodies [14-16] has been conducted in the VST as well. 

Figure 2 depicts a diagram of this VST. The VST is an atmospheric closed-loop, annular 
return wind tunnel. It features a 12-sided test section that is 20 feet in diameter and 25 feet 
tall. The test section velocity can be varied from 0 to approximately 85 ft/s. The airflow is 
produced by a 400 hp DC electric motor with a 3-bladed fixed pitch fan. A tether can be 
attached to the test article to prevent wall impacts, set initial conditions, and provide quick 
recovery. 

The Tunnel Coordinate System (TCS) depicted in Figure 2 has been adopted as the 
inertial reference frame for analyzing the test data. The origin is located in the center of the 
test section. The X-axis is aligned with the central axis of the tunnel pointed toward the 
floor. The Y-axis and Z-axis are pointed in arbitrary directions with both axes orthogonal 
with the X-axis and each other. 



Figure 2. Cross-sectional diagram of the 20- ft Vertical Spin Tunnel. 
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2.2 Test Articles 


Two different model configurations were flown inside the VST during this experiment. 
The first configuration has the canards on the LAV (see Figure 1) in a stowed state. The 
second configuration has the canards deployed. 

Table 1 provides a summary of the model prop- Table 1. Summary of physical prop- 
erties used for this analysis. The test models used in erties for the test models. 


the investigation are 1/29-th dynamically scaled rep- 


resentations of the current LAV design at the time of 

Property 

Value 

the experiment. The dynamic scaling was performed 

L , m 


0.173 

using the method outlined by Wolowicz, Bowman, 

S, m 2 


0.0236 

and Gilbert [17] for free flight test articles in the in- 

m , kg 


0.6187 

compressible flow regime. Target values for the test 

I XX 1 

kg- 

2 

-m 

0.00100 

vehicles mass, Moments of Inertia (MOI), and Prod- 

lyyi 

kg- 

2 

■m 

0.01488 

ucts of Inertia (POI) are set based on scaling the full 

I ZZ 5 

kg- 

2 

m 

0.01471 

size LAV. The test vehicles are ballasted to within 

Ixy 5 

kg- 

2 

■m 

0.00009 

±2% of the target values. 

lyzi 

kg- 

2 

m 

0 

Figure 3 describes the test article’s Body Coordi- 

1-XZ 5 

kg- 

2 

■m 

-0.00018 

nate System (BCS) (Xb, Yb, Zb) with the Center 





of Mass (CM) at the origin. The moment coefficients (Cj, 

Cm i 

c n 

) described in Figure 


3 are standard. For this analysis, the aerodynamc force acting on the LAV is described 
by the force coefficients based on the LAVcoordinate system {Cx, Cy , Cz)- This system 
was chosen for its simplicity regarding the data reduction process, a describes the angle of 
attack. /3 is the side-slip angle. 



Figure 3. Body Coordinate System for the Launch Abort Vehicle. 
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2.3 Data Acquisition System 

The Data Acquisition System (DAS) for this experiment consisted of three parts. The 
first part is a customized commercially available motion capture system for tracking the 
test article. The second part is a pitot-static system measuring the dynamic pressure in the 
wind tunnel. The third part is a commercial video camera for visual recording of the test. 
The data presented in this report comes from the motion capture and pitot-static systems. 

The motion capture system utilizes an array of 
8 cameras to capture images of the test section from 
different vantage points. Small reflective targets, Fig- 
ure 4 [18], on the model are used as reference targets 
for the DAS to compare image frames. Using pho- 
togrammetry [19,20], the position and orientation of 
the model relative to the TCS is obtained. A more 
detailed description of photogrammetry can be found 
in Reference [20, www.geodetic.com]. 

The frame rate of the motion capture in these 
experiments is 150 Hz. The data returned by the 
system are the CM position (X, Y , Z ) and Euler 
angles ( <f > , 0, t/j). The Euler angles correspond to a 
X-Y-Z rotation sequence to map TCS measurements 
to the BCS. 

A pitot-static tube measures dynamic pressure 
( q ) in the wind tunnel. The VST is atmospheric and operates in the subsonic incompress- 
ible flow regime. A barometer and thermometer located at the VST are used to measure 
ambient air pressure and temperature so that air density may be calculated. The air density 
is assumed to be constant through the VST during the testing. The tunnel airspeed (V) is 
calculate using the measured dynamic pressure and calculated air density. 

It should be noted that the motion capture system has been recently updated using 
equipment developed by Vicon®. A rigorous error analysis has not been conducted on the 
system. An informal check regarding the precision of the system is in Reference [18]. The 
results obtained from this check can be seen in Table 2. There was no attempt to quantify 
the uncertainty of the wind tunnel velocity and dynamic pressure measurements. 



Figure 4. Reflective targets used to 
capture model motion. 


Table 2. Precision of data measured by the Vicon® data acquisition system. 


Measurment 

a 

<f>, rad 

0.0002 

9, rad 

0.0002 

V’, rad 

0.0003 

V, m 

0.0001 

Y , m 

0.0001 

Z , m 

0.0001 
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3 Analysis Method 


This section discusses the approach used to model the aerodynamic forces and moments 
acting on the LAV. Section 3.1 describes the data reduction process for obtaining useful 
aerodynamics information from experiment observations. Section 3.2 discusses the methods 
used to obtain the model structures and parameter estimates for the aerodynamic forces 
and moments. 


3.1 Data Reduction 

SIDPAC was designed for aircraft which 
could be fully instrumented with sensors 
(rate gyros, accelerometers, etc) to collect 
data. However, the data collected for this 
investigation was limited to data concerning 
the model’s CM position (X, Y, Z), orienta- 
tion (Euler angles <f>, 9, ip), dynamic pressure 
(q), and airspeed( V) inside the VST. The 
first task of this analysis is to extract the 
aerodynamic state, forces, and moments. 
Figure 5 outlines the process used to extract 
this information from the experiment ob- 
servations. Subsections 3.1.1 though 3.1.7 
details the methods used to obtain the in- 
formation described in Figure 5. 

3.1.1 Data Smoothing 



Figure 5. Schematic of data reduction process 
used to extract aerodynamic information from 
the experiment observations. 


Data smoothing removes low power high frequency content from the data. The high 
frequency content usually contains either noise or higher-order structural responses which 
are outside the focus of this research. Smoothing data prior to differentiation prevents large 
increases in noise [21]. 


The smoothing scheme used for this re- 
search uses Fourier analysis to treat the 
smoothing problem in a global sense. This 
method is outlined by Morelli and Klein 
in Reference [22, Section 11.2.3 Global 
Smoothing in the Frequency Domain ]. This 
method requires the analyst to select cut- 
off frequencies to smooth the signals. These 
cut-off frequencies need to be large enough 
to capture the low frequency dynamics and 
small enough so that high frequency noise 
is removed. The cut-off frequencies used 
in this analysis were placed where the fre- 
quency content began to decrease towards 
zero. Figure 6 indicates where cut-off fre- 
quencies are chosen relative to the fre- 
quency content observed in the VST data. 
The SIDPAC algorithm, “smoo” [5, 22] is 
used to implement this smoothing technique. 



Figure 6. Cut-off frequency selection for 
smoothing. 


6 


3.1.2 Transformation Matrices 


Two transformation matrices are required for analyzing the data. One transformation 
matrix is the Direction Cosine Matrix (DCM), [C], to transform vectors from the TCS, see 
Figure 2 on page 3, to the BCS, see Figure 3 on page 4. The kinematic transformation 
matrix, [T], is used to transform the Euler rates into the angular velocity vector in the body 
frame. 

The DCM is constructed using the Euler angles for the X-Y-Z rotation sequence [23-26]. 
The first rotation is for roll, </>. The second rotation is pitch, 9. The third rotation is yaw, 
ip. Eq. 1 (c is for cosine, s is for sine) yields the DCM that transforms data from TCS to 
BCS at each point in time. 


[C] 


cQap cips9s<p + sipccp 
—c9sip —sips9s<p + cipccp 
s9 — c9s(j> 


—cips9c(p + sipsep 
sips9c<p + cipscp 
c9c<p 


(1) 


[T] is required to transform the Euler angle rates into body angular rates. Eq. 2 is 
derived from the X-Y-Z DCM rotation sequence. 


[T\ 


cos 9 cos ip 
— cos 9 sin ip 
sin0 


sin^> 

cosip 


0 


0 

0 

1 


( 2 ) 


3.1.3 Inertial and Body Velocities 

The inertial velocity vector is calculated by differentiating the position data with respect 
to time. A five point running 2 nd order polynomial performs this differentiation. This 
method calculates the derivative at each time step using the current data point and its 4 
closest neighbors. The four closest neighbors includes two past and two future data points 
for the majority of the data. The 4 neighboring points become one-sided (either all past 
or all future data points) as the end points of the data stream are approached. Data close 
to the endpoints are typically left out of the analysis. A more detailed description of this 
method can be found in Reference [22, Section 11.4 Numerical Differentiation \. 

Once the inertial velocity is determined, the body velocity can be calculated using the 
DCM matrix for each data point. Eq. 3 describes the transformation. Angular velocity 
does not appear in this transform because the CM is collocated with the BCS origin (refer 
to Figure 3 on page 4). 


u 


' x ' 

V 

= \c\ 

Y 

w 


z _ 


3.1.4 Euler and Body Angular Rates 


The Euler angles are differentiated using the running polynomial discussed in Section 
3.1.3 to determine the Euler rates (<p, 9, ip)- Eq. 4 transforms the Euler rates into body 
angular velocity components. 


p 


’ (P ’ 

q 

= [T] 

9 

r 


. ip . 
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3.1.5 Aerodynamic Orientation 


The aerodynamic orientation describes the direction of the relative wind vector. The 
relative wind vector (u 1 ,v' ,w') (See Eq. 5) is a combination of the vehicle’s velocity (u, v , 
w) and the tunnel’s free stream vector in the BCS. The tunnel’s free stream is assumed to 
be acting uniformly in the TCS negative x-direction with magnitude V (refer to Figure 2 
on page 3). 


' v! ' 


u 


' -V ' 

v' 

= 

V 

-[ci 

0 

w' 


w 


0 


Vehicle Tunnel 


( 5 ) 


The vehicle’s velocity is small in comparison to the wind tunnel free stream vector and can 
be neglected. The calculation for the relative wind vector components are simplified to Eq. 
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v! = Vcos9cosil’ (6a) 

v 1 = —V cos 9 sin ip (6b) 

w' = V sin 6 (6c) 

The relative wind vector’s orientation is typically given in terms of incidence angles. 
Angle of attack and sideslip angle are typically used for this purpose. The angle of attack 
and sideslip angle can be calculated directly using Eqs. 7 and 8. 

w' 

tana = — (7) 

u 

sin/3 = V - (8) 

Total angle of attack and the aerodynamic clock angle are a different set of incidence 
angles which describe the relative wind vector’s orientation. The total angle of attack, olt, 
is the angle between the relative wind vector and the body x-axis component of the wind 
vector (Eq. 9). The aerodynamic clock angle, <pA , is the angle between the body y-axis and 
z-axis components of the relative wind vector (Eq. 10). 

v! , . 

cosar = — (9) 


tan 4>a = (10) 

w' 

NOTE: Most arctangent algorithms provide solutions between —tt/2 and n/2 (Cartesian 
quadrants 1 and 4) because they cannot distinguish quadrants 2 and 3. It is important 
to use component based arctangent algorithms to overcome this deficiency. This research 
uses an algorithm that provides solutions between 0 and 2-7 t for angle of attack. A different 
algorithm provides solutions between — ir and 7 r for the aerodynamic clock angle. 


3.1.6 Inertial Accelerations and Aerodynamic Forces 

The inertial acceleration of the test article (A, Y , Z ) are obtained by taking the time 
derivative of the inertial velocity. The time derivatives were taken using the same method 
described earlier for the inertial velocity. This information can then be used to determine 
the aerodynamic forces acting on the body. The aerodynamic forces acting on the model 
are calculated using the Eq. 11 (c is for cosine, s is for sine). For this research the forces 
are presented in terms of the coefficients (Cx, Cy, Cz) acting along the BCS axes (refer to 
Figure 3 on page 4). 


Cx = ^ 
qS 

— g^j cdcijj + Y (ci/isdsc/) + sipccf) + Z {—cijisdccj) + sipscj)) 

s|% 

II 

& 

— ^ X — g'j c9si[) + Y (—sipsdscj) + c?/>c^>) + Z ( stpsOccf) + cipscj)] 

£ I^S 
ii 
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3.1.7 Body Angular Acceleration and Aerodynamic Moments 

The body angular accelerations are obtained by differentiating the body angular velocities 
using the method presented in Section 3.1.3. Euler’s equations of motion for rigid bodies can 
now be applied to the model. These equations describe the net moment acting on the body 
in the BCS. The moments can be calculated using the body angular rates and accelerations. 
Eq. 12 displays Euler’s equations in terms of the aerodynamic moment coefficients (Cj, C mi 
C n ) acting on the model. 


Ci 


C n 


C n 


\jxxP IxyQ l'x z f ’ “t“ lyy) “b lyzi?* Q ) T IxyP'T’ I xzQp\ 

\^yyQ ^xyP lyzf T P^ijxx Izz) T Ixz{p C ) T lyzQP IxyQf\ 

-j- [izzf - IxzP - Iy Z q + qp{Iyy - hx) + Ixy (<? 2 - p 2 ) + I XZ qr - I yz pr \ 


(12a) 

(12b) 

(12c) 


3.2 Data Analysis 

This section discusses the tools used to analyze the data and determine models for the 
aerodynamic forces and moments. Stepwise regression is performed in both the time and 
frequency domains to determine the model structure for the forces and moments. The 
method of least squares is used to determine parameter estimates and uncertainties for the 
force and moment models for each test case. Mean models for each force and moment are 
developed using covariance weighted averages of the parameter estimates and covariances 
found for individual test flights. 

Sections 3.2.1 through 3.2.4 discusses how these tools are used in more detail. Section 
3.2.1 describes the least square method and implementation. Section 3.2.2 describes the 
tools used to transform the variable time histories in to the frequency domain. Section 3.2.3 
describes the stepwise regression process and how it is implemented in this research. Section 
3.2.4 describes the covariance weighted mean. 
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3.2.1 Ordinary Least Squares 


Ordinary least squares analysis is the method used to obtain parameter estimates and 
uncertainties. The stepwise regression uses the least squares model shown in Eq. 13. H is 
a matrix of regressors and x is a vector of unknown parameters, v is the error between the 
model and measurements. Eq. 14 describes the estimate of x using the observations and 
regressors matrix. H* is the conjugate transpose of H. 

z = Hx + v (13) 

x = [Re (H*H)] -1 Re (H*z) (14) 

Ordinary least squares provides a covariance matrix on the estimates. This covariance 
can be used to determine the variance on the individual parameters and the correlation 
between the parameters. The covariance matrix is determined using Eq. 15. 

r * = ( 15 ) 


To hasten the analysis process, the full data streams for each dependent variable and 
regressor are not used. Short segments away from the endpoints are used to analyze the data. 
During flight tests, model residuals are often correlated with neighboring data points [22], 
These correlations have a serious effect on the covariance matrix calculated for parameter 
estimates. Reference [22, Section 5.2, Generalized Least Squares] provides an in-depth 
discussion on the post processing needed to correct the data. The SIDPAC algorithm 
“r_colores.m” [5,22] is used to compute the parameter estimates and covariance matrix for 
model parameters. This algorithm account for the frequency content of the residuals as well 
as providing a least squares fit to the data. 

3.2.2 Frequency Domain 

Performing the analysis in the frequency domain as well as the time domain is a useful 
analysis technique. Transforming the dependent variable and regressors into the frequency 
domain helps validate model parameters. Model parameters need to be consistent between 
both the time and frequency domains. This is due to the independent nature of each 
regressor. 

The z-transform [22] performs the Fourier transform with improved computational effi- 
ciency and arbitrary frequency selection. This method uses interpolation of data points to 
approximate a function with respect to time. This decouples the available frequencies from 
the data points. The Discrete Fourier Transform (DFT) is then utilized only on the desired 
band of frequencies. 

This feature is very helpful if the data is oversampled. The data presented in this research 
is sampled at 150 Hz. The standard DFT would compute the frequency domain coefficients 
for each frequency up to the Nyquist frequency of 75 Hz. The main signal content presented 
in the data is typically below 5 Hz (See Figure 6). 

3.2.3 Stepwise Regression 

Stepwise regression will be applied to the force and moment coefficients obtained from 
the data reduction discussed previously. Stepwise regression is a method for identifying the 
model structure of a regression equation using a pool of candidate regressors (also obtained 
from the data reduction) . The model structure represents the number and form of terms in 
the regression equation. A more thorough discussion of stepwise regression can be found in 
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Reference [22, Section 5.4.2 Stepwise Regression], The SIDPAC algorithm “swr.m” is used 
to perform the regression. This algorithm uses the method of least squares to determine 
parameter estimates for the model. 

In this research, pools of candidate regressors were formed for each force and moment 
coefficient. Models are initially formed with only a constant term. A step is taken when a 
single regressor is either added or removed from the model. Parameter estimates for each 
regressor and the covariance matrix of the model are updated after each step. The F-ratio 
(see Eq. 16) is calculated for each included regressor to measure the regressors’ importance. 
Correlations between model residual and excluded regressors are recomputed. 



(16) 


Regressors are added to the model beginning with the largest correlation until either all 
regressors are added or the largest correlation with the remaining excluded regressors is less 
than 0.01. The F-ratio is then checked to see if model parameters need to be removed from 
the model. A minimum F-ratio of 20 (5 times the value for 95% confidence) is required for 
all regressors used in modeling the data for this investigation. Regressors that satisfy the 
minimum F-ratio can also be excluded depending on how the regressor’s F-ratio compares 
to the F-ratio of other regressors and how much the coefficient of determination (R 2 ) of 
the model changes when the regressor is added or removed from the model. Regressors 
with F-ratios substantially lower than other regressors’ F-ratios will not contribute much 
information to the model and can be neglected. This is seen by tracking how much a 
model’s R 2 changes when a regressor is added or removed. If the change in R 2 is small 
(< 2% ) and F-ratio is small compared to others, then the regressor may be removed 
without degrading the model. A more complete discussion of the F-ratio and it’s uses can 
be found in Reference [22, Section 5.4.3 Stepwise Regression]. 


3.2.4 Covariance Weighted Mean 

This research works with multiple data segments from different test runs. While the 
model structure may be similar between test runs, the parameter estimates for each test 
will slightly vary between test runs. A covariance weighted mean is used to combine the 
different sets of parameter estimates into a single average set to describe all the test data. 

Eqs 17 and 18 [27] describes the mean estimate and its covariance, respectively, for n 
tests. x ? ; and T,; represent the parameter estimates and covariance matrix for each test. 


x = 



r, = 



(17) 

(18) 


Using the covariance matrix as the weighting factor places emphasis on the parameter 
estimates with low uncertainty. Care must be taken with covariance matrices that have 
large correlation terms between parameter estimates. Large correlations can move the mean 
parameter estimates outside of the values provided by the separate estimates. 
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4 Results and Discussion 


This section contains the results of the experiment data analysis. Models were identified 
for C m , C n , Cy , and Cz- A detailed discussion of these models can be found in sections 
4. 1-4.4. A summary of the models and parameter estimates for these force and moment 
coefficients are provided in Section 4.5. 

4.1 Pitch Moment Coefficient 

Figure 7 shows a plot of C m versus a for the LAV in the stowed-canard configuration for 
tests 1, 2, and 3. Static wind tunnel data obtained by the TDT (Ref. [28]) is compared to 
the VST data in Figure 7. It must be noted that the TDT data was obtained at a Mach 
number of 0.2 and Reynold’s number of 4.7e+6. The VST data was obtained at Mach 
0.06 and Reynold’s number of 3.0e+5. An “apples-to-apples” comparison cannot be made 
between the two datasets. However, the comparison is still useful for spotting trends in the 
data. The final model for C m is also shown in Figure 7 for comparison. 



'140 150 160 170 180 190 200 210 220 
a, deg 



140 150 160 170 180 190 200 210 220 
a, deg 



140 150 160 170 180 190 200 210 220 
a, deg 


(a) Test 1. 


(b) Test 2. 


(c) Test 3. 


Figure 7. C m versus a comparison between Vertical Spin Tunnel (VST) data, Transonic 
Dynamics Tunnel (TDT) data, and the C m model for the stowed-canard configuration. VST 
data is at M = 0.06 and Re = 3e+5. TDT data is at M = 0.2 and Re = 4.7e+6. Error 
bars represent ler standard deviation. 


Figure 7 indicates that the LAV is statically stable in the heat shield-forward flight 
direction. Figure 7(a) indicates a highly nonlinear relationship between a and C m for test 
1, while Figures 7(b) and 7(c) demonstrate a more linear relationship with some nonlinear 
trend observed at the extremes of the a range. The VST data and model shown in Figure 
7 describe a steeper slope than the TDT data. Figure 7(b) and 7(c) indicate that VST and 
TDT data are in close agreement at angles of attack near 180 degrees. However, the VST 
and TDT diverge at angles of attack further away from 180 degrees. The C m — a slope 
appear to lessen toward the extremes of the TDT data. The VST data indicate that the 
C m — a slope increase near the extremes. 

Figure 8 describes the C m -a behavior for the LAV in the deployed-canard configuration 
seen in tests 4 and 8. The identified model is displayed with the VST data for comparison. 
No static wind tunnel data was available for comparison for the deployed-canard configura- 
tion. Test 4 in Figure 8(a) shows a slightly cubic relationship between C m and a. However, 
the data also appears to lie along the line representing the identified model for the deployed- 
canard configuration. Figure 8(b) (test 8) demonstrates a more linear relationship between 
Cm and a which almost falls on the model line exactly with a couple outliers. Figure 8 shows 
that the zero crossing is shifted from 180 degrees to about 174 degrees angle of attack. 
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(a) Test 4. 


(b) Test 8. 


Figure 8. C m vs a comparison between Vertical Spin Tunnel (VST) data and the C m model 
for the deployed-canard configuration. Vertical Spin Tunnel (VST) data is at M = 0.06 and 
Re = 3e+5. Error bars represent la standard deviation. 


Table 3. C m model structures identified using stepwise regression. 


Tests Time Domain Frequency Domain 
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“Stowed-canard configuration. 
b Deployed-canard configuration. 


Table 4. Parameter estimates for C m . 



Time Domain 

Frequency Domain 

Tests 

Cm 0 i 0 ’) 

Cm a , (<r) 

C ma , (a) 

l a 

-0.007 (0.002) 

-0.50 (0.01) 

-0.51 (0.01) 

2 a 

-0.005 (0.002) 

-0.50 (0.01) 

-0.497 (0.008) 

3 a 

-0.010 (0.004) 

-0.53 (0.01) 

-0.53 (0.02) 

Mean 0 

-0.006 (0.001) 

-0.513 (0.007) 

-0.505 (0.006) 

4 b 

-0.076 (0.003) 

-0.67 (0.01) 

-0.67 (0.01) c 

8 b 

-0.063 (0.003) 

-0.615 (0.009) 

-0.62 (0.02) c 

Mean 6 

-0.069 (0.002) 

-0.631 (0.008) 

-0.65 (0.01) 


a Stowed-canards configuration. 
^Deployed-canards configuration. 
c Frequency domain estimate adjusted for bias. 
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Table 3 provides a summary of the model structures identified by the stepwise regression 
in the time and frequency domains. All of the identified models feature a Cm a , term for 
both configurations. C m , represents the slope of C m — a near the angle of attack at 180 
degrees ( 7 r radians). A regressor of (a — n) was chosen to represent the data since the 
angle of attack data stayed close to 180 degrees. Only the regressions of tests 1 and 4 in 
the time domain indicate that the cubic angle of attack term is significant for the canard- 
stowed and deployed-canard configurations, respectively. This makes sense given the cubic 
relationship between C m — a seen in Figures 7(a) and 8(a), respectively. This cubic term 
may be representative of wind tunnel effects that are not captured by the data. These 
unobserved effects may be responsible for why test 1 remained in the same Y-Z location for 
the duration of its flight unlike tests 2 and 3 which moved about the VST. 

Table 3 shows that the stepwise regression indicates that the body angular roll rate ( p ) is 
significant to the deployed-canard configuration model in the frequency domain. This is due 
to the increased constant term ( C mo ) when the canards are deployed. This was proven by 
calculating C m — C mo in the frequency domain and redoing the stepwise regression for tests 
4 and 8. Removing the constant term in the time domain decreased the importance of p in 
the frequency domain. The final model for C m is shown in Equation 19. The terms shown 
in Equation 19 are consistent throughout the tests. Equation 19 is simple and accounts for 
most of the data obtained through the tests. 

Cm = Cm 0 + Cm a i (<T _ ?r) (19) 

The stepwise regression reveals Equation 19 to be the common model between all tests 
described in Table 3 for both stowed-canard and deployed-canard configurations. Table 4 
shows the parameter estimates and the covariance weighted mean estimates for the two 
configurations in the time and frequency domains. Figure 18 provides a graphical compar- 
ison of the parameter estimates described in Table 4. Figure 9 shows that both time and 
frequency domain estimates for C ma , are nearly identical. The estimates for C mo and C ma , 
do not significantly change between test runs which indicates an accurate model structure 
was identified. 
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(b) Cm, 


Figure 9. Comparison of pitch moment coefficient model parameters. Error bars represent 
ler standard deviation. 
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4.2 Yaw Moment Coefficient 


Figure 10 shows a plot of C n versus j3 for the stowed-canard configuration for tests 1, 2, 
and 3. Figure 10 compares the VST data to the identified C n model. Figure 10 indicates 
that the LAV is statically stable in the heat shield-forward flight direction. Figure 10(a) 
shows a very linear relationship between j3 and C n for test 1 that falls almost exactly on 
the model line. Figures 10(b) and 10(c) indicate a slightly nonlinear relationship between 
C n and /3 for tests 2 and 3. The linear C n model still appears to accurately reflect the VST 
data captured during testing. 



-40 -30 -20 -10 0 10 20 30 40 

P, deg 


(a) Test 1. 



-Uo -30 -20 -10 0 10 20 30 40 

P, deg 

(b) Test 2. 



-40 -30 -20 -10 0 10 20 30 40 

P, deg 

(c) Test 3. 


Figure 10. C n versus /3 comparison between Vertical Spin Tunnel (VST) data and the C n 
model for the stowed-canard configuration. VST data is at M = 0.06 and Re = 3e+5. Error 
bars represent la standard deviation. 

Figure 11 describes the C n -j3 behavior for the deployed-canard configuration seen in 
tests 4 and 8. The identified model is displayed with the VST data for comparison. No 
static wind tunnel data was available for comparison for the deployed-canard configuration. 
Figure 10 shows that both tests show a linear between C n and /3 which is captured by the 
identified model. However, the data also appears to lie along the identified model for the 
deployed-canard configuration. 



Uo -30 -20 -10 0 10 20 30 40 

P, deg 

(a) Test 4. 
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P, deg 

(b) Test 8. 


Figure 11. C n versus /3 comparison between Vertical Spin Tunnel (VST) data and the C n 
model for the deployed-canard configuration. VST data is at M = 0.06 and Re = 3e+5. 
Error bars represent ler standard deviation. 

Table 5 provides a summary of the model structures identified using stepwise regression 
in the time and frequency domains. The time and frequency domain model structures consis- 
tently describe a linear model (Equation 20). C no is the constant bias term and only appears 
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Table 5. C n model structures identified using stepwise regression. 


Tests 

Time Domain 

Frequency Domain 

l a 

C no + C n/3 f3 

CnpP 

2 a 

Cno + C n , f3 

C np f3 

3“ 

Cno "F C n , (3 

Cnp (3 

4 b 

Cn 0 "F C n g (3 

CnpP 

8 b 

Cno "F C n/3 f3 

Cnpf3 


“Stowed-canard configuration. 
^Deployed-canard configuration. 


in the time domain. Table 6 describes the least squares parameter estimates and covariance 
weighted mean estimates for both stowed-canard and deployed-canard configurations. 

C n = C no + C„,/3 (20) 

Figure 12 provides a graphical comparison of the parameter estimates for C no and C nfj . 
The time domain estimates for C no have nearly the same value for all the stowed canard 
test runs. C„ 0 appears to significantly change between test runs for the deployed canard 
configuration. The reason the deployed canard mean value is closer to the test 8 is because 
of the smaller uncertainty and a significant correlation between C n/3 and C no which would 
drive the covariance weighted mean toward the test 8 estimate. This same phenomena can 
be seen in C np for the stowed canard test runs. 


Table 6. Parameter estimates for C n . 



Time Domain 

Frequency Domain 

Tests 

Cn 0 (*) 

Cn, (a) 

Cn, (a) 

1“ 

0.008 (0.002) 

-0.59 (0.01) 

-0.59 (0.01) 

2“ 

0.009 (0.001) 

-0.51 (0.02) 

-0.507 (0.007) 

3“ 

0.013 (0.003) 

-0.54 (0.01) 

-0.54 (0.01) 

Mean a 

0.0092 (0.0008) 

-0.561 (0.007) 

-0.532 (0.005) 

4 b 

0.008 (0.002) 

-0.62 (0.01) 

-0.61 (0.02) 

8 b 

-0.007 (0.001) 

-0.56 (0.01) 

-0.56 (0.02) 

Mean 6 

-0.0057 (0.0009) 

-0.586 (0.008) 

-0.59 (0.01) 


“Stowed-canard configuration. 
b Deployed-canard configuration. 
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(b) C„ 


Figure 12. Comparison of yaw moment coefficient model parameters. Error bars represent 
la standard deviation. 


4.3 Body Y-Force Coefficient 

Figure 13 shows a plot of Cy versus /3 for the LAV in the stowed-canard configuration 
for tests 1, 2, and 3. Figure 13 displays the model (Equation 21) for Cy without the Cy r 
term. Static wind tunnel data was not explicitly available for comparison. The VST data is 
displayed based on the average Cy for different partitions of side slip angle (/?). The error 
bar represents ler uncertainty based on the partition populations. 



-40 -30 -20 -10 0 10 20 30 40 

P, deg 
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(a) Test 1. 


(b) Test 2. 



P, deg 


(c) Test 3. 


Figure 13. Cy versus /? comparison between Vertical Spin Tunnel (VST) data and the Cy 
model for the stowed-canard configuration. VST data is at M = 0.06 and Re = 3e+5. Error 
bars represent ler standard deviation. 

Figure 13(a) indicates that Cy has a very linear relationship with /3. Figures 13(b) and 
13(c) describe a slightly nonlinear relationship between Cy and /3. The linear model appears 
to track Cy very well and start to diverge around the endpoints. 

Figure 14 describes the Cy-/3 behavior for the LAV in the deployed-canard configuration 
seen in tests 4 and 8. The model is calculated using Equation 21 without the Cy r term. 
Test 4 in Figure 14(a) describes a linear Cy — /3 relationship for test 4. Test 8 (Figure 14(b)) 
shows a nonlinear relationship between Cy and /3. 

Table 7 provides a summary of the model structures identified by the stepwise regression 
in the time and frequency domains. Test 3 is the only test in which the stepwise regression 
yields a model structure that is different from Equation 21. Cy p accounts for the majority 
of the frequency content seen in the data. While Cy r is identified by the stepwise regression, 
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(b) Test 8. 


Figure 14. Cy versus (3 comparison between Vertical Spin Tunnel (VST) data and the Cy 
model for the deployed-canard configuration. VST data is at M = 0.06 and Re = 3e+5. 
Error bars represent ler standard deviation. 


it tends to be a much smaller contributer to the model. 


Table 7. C n model structures identified using stepwise regression. 


Tests 

Time Domain 

Frequency Domain 

1“ 

Cyo+CypP + Cy^ 

c Yf) P + c Yr £ 

2“ 

c Yo + c Yp p + c Yr £ 

c Yp p + c Yr $ 

3“ 

Cy 0 +Cy^ + C Yq § 

Cr^ + C Yq § 

4 b 

c Yo + c Yp p + c Yr ^ 

C Yp l3 + C Yr £ 

8 b 

c Yo + c Yp p + c Yr $ 

c Yp p + c Yr £ 


“Stowed-canard configuration. 
^Deployed-canard configuration. 


Table 8 displays the parameter estimates and the covariance weighted mean estimates 
for the two configurations in the time and frequency domains. The values are the results 
of the ordinary least squares method described in Section 3.2.1. Notice that the frequency 
domain values for Cy r are not as consistent as the time domain values. This may indicate 
that only C Yp is necessary for an adequate model. 


Cy = Cy 0 + C Y J + C Yr — (21) 

Table 8 indicates that the canards have a minor affect on Cy. Cy B and Cy r appear to 
change only slightly based on the mean estimates. There is a slight increase in the magnitude 
of C Yp , and a small decrease in the magnitude of Cy r . 

Figure 15 displays a graphical comparison of the model parameters in Table 8 for each 
test and mean combination for the stowed and deployed canard configurations. Figure 
15(a) indicates that C yo may be sensitive to the test run, with the estimate for test 3 being 
significantly different from the other tests. Figures 15(b) and 15(c) display frequency domain 
estimates with large uncertainties on parameters C Yp and Cy T for test 3. This may signal 
that the identified model may not be sufficient for test 3 and that some information was not 
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Table 8. Parameter estimates for Gy 



Time Domain 


Frequency Domain 

Tests 

Cy 0 (a) 

Cy„ W 

C Yr (a) 

Cy p (*) 

C Yr (a) 

1“ 

0.001 (0.003) 

-0.468 (0.009) 

-2.5 (0.3) 

-0.48 (0.08) 

-4 (2) 

2“ 

0.004 (0.001) 

-0.396 (0.008) 

-1.9 (0.2) 

-0.34 (0.04) 

-2(1) 

3“ 

-0.014 (0.003) 

-0.38 (0.02) 

-2.0 (0.5) 

-0.5 (0.1) 

-3 (4) 

Mean 6 

0.001 (0.001) 

-0.422 (0.006) 

-2.2 (0.2) 

-0.37 (0.03) 

-3.5 (0.9) 

4 6 

-0.0017 (0.0008) 

-0.493 (0.007) 

-1.4 (0.2) 

-0.50 (0.04) 

-2(1) 

8 b 

-0.001 (0.003) 

-0.40 (0.03) 

-2.0 (0.7) 

-0.35 (0.07) 

-2 (2) 

Mean 6 

-0.0016 (0.0007) 

-0.487 (0.007) 

-1.5 (0.2) 

-0.46 (0.03) 

-1.5 (0.9) 


“Stowed-canard configuration. 
^Deployed-canard configuration. 


captured by the data. The deployed canard configuration estimates remain similar though 
out the test runs indicating an accurate model was identified. 
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Figure 15. Comparison of Y-force coefficient model parameters. Error bars represent lcr 
standard deviation. 
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4.4 Body Z-Force Coefficient 

Figure 16 shows a plot of Cz versus a for the LAV in the stowed-canard configuration 
for tests 1, 2, and 3. Static wind tunnel data obtained by the TDT is compared to the VST 
data in Figure 16. The identified model is also plotted alongside the TDT and VST data for 
comparison. The TDT data was obtained at a Mach number of 0.2 and Reynold’s number 
of 4.7e+6. The VST data was obtained at Mach 0.06 and Reynold’s number of 3.0e+5. The 
model uses only Cz 0 and Cz , terms. 
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(a) Test 1. 


(b) Test 2. 


(c) Test 3. 


Figure 16. Cz versus a comparison between Vertical Spin Tunnel (VST) data, Transonic 
Dynamics Tunnel (TDT) data, and the Cz model for the stowed-canard configuration. VST 
data is at M = 0.06 and Re = 3e+5. TDT data is at M = 0.2 and Re = 4.7e+6.Error bars 
represent ler standard deviation. 


Figure 16 showes that the TDT and VST present different Cz — a behaviors at angles 
of attack close to 180 degrees. The TDT data presents a highly nonlinear relationship 
compared to the VST data. This difference is most likely due to the difference in Reynolds 
numbers. It is interesting to note that the slope of the Cz versus a for the TDT data 
becomes similar to the slope represented by the VST data at angles of attack farther away 
from 180 degrees. 
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(a) Test 4. 


(b) Test 8. 


Figure 17. Cz versus a comparison between Vertical Spin Tunnel (VST) data and the Cz 
model for the deployed-canard configuration. VST data is at M = 0.06 and Re = 3e+5. 
Error bars represent ler standard deviation. 


Figure 17 describes the Cz-a behavior for the LAV in the deployed-canard configuration 
seen in tests 4 and 8. The identified model displayed with the VST data for comparison. No 
static wind tunnel data was available for comparison for the deployed-canard configuration. 
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Test 8 in Figure 17(b) shows a slightly nonlinear relationship between Cz and a. Figure 
17 shows that the zero crossing is shifted from 180 degrees to about 177 degrees angle of 
attack. 

Table 9 provides a summary of the model structures identified by the stepwise regression 
in the time and frequency domains. All of the identified models feature the model (Equation 
22) for both configurations. Only test 8 appears to favor an additional nonlinear term 
(Cz l3 )- This makes sense given the cubic relationship between Cz — a seen in Figure 
17(b). No other test indicates a strong cubic relationship between Cz and a. 


Table 9. Cz model structures identified using stepwise regression. 


Tests 

Time Domain 


Frequency Domain 

l a 

Cz 0 + C z 

( a ~ 7r ) + Cz q fp- 


Cz a , 

(a-n) + C Zg £ 

2° 

Cz 0 + Cz t 

„/ ( a — 7r) + C Zq 


C Za , 

(a — 7r) + C Zq 

3“ 

Cz 0 + Cz a , (a 

- 7r) + C Zq fp- + Cz r Tff 

Cz a , 

(a- 

tt) + Cz q fp- + Cz r 

4 b 

Cz 0 + C z , 

„/ ( a ~ 7r ) + Cz q fp- 

C Za , 

(a - 

tt) + C Zp fp- + C Zq 

8 b 

C Zo + Cz a , (a — tt) + Cz a , 3 (a - tt) 3 

C z , 

(«- 

- tt) + C Za , 3 (a -n) 3 



+C Zq w 


+Cz p §£ + C Zq fp- 


a Stowed-canard Configuration. 
^Deployed-canard configuration. 


Table 9 shows that the stepwise regression indicates that the body angular roll rate (p) 
is significant to the deployed-canard configuration model in the frequency domain. This 
is likely due to the same reasons why p appeared in the frequency domain regression for 
C m . Equation 22 also displays a symmetry with Equation 21 which makes sense given the 
axisymmetric nature of the LAV. 


Cz = C Zo + Cz a , (a — n) + C Zq y ( 22 ) 

Table 10 shows the parameter estimates that were obtained using Equation 22 for the 
model structure. Values are rounded based on the significant figures indicated by their 
uncertainty (a). The canards appear to increase the values of the model parameters. Cz 0 
increases by about 0.02, Cz a , increases by 0.1, and Cz q increases by about 0.5. 

Figure 18 displays a graphical comparison of the model parameters in Table 10 for each 
test and mean combination for the stowed and deployed canard configurations. Figure 18(a) 
shows consistent estimates for Cz 0 for both the stowed and deployed canard test runs. Figure 
18(b) shows similar consistent time and frequency domain for test runs with the exception of 
test 3. Test 3 provided a much lower estimate for Cz a , with larger uncertainty. Figure 18(c) 
shows consistent estimates for Cz q for both stowed and deployed canard configurations. 
However, the uncertainties on Cz q for tests 1 and 3 indicate that some information may be 
missing from either the model structure or test data. 
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Table 10. Parameter estimates for Cz- 



Time Domain 


Frequency Domain 

Tests 

Cz 0 (a) 

Cz a , (a) 

Cz, (a) 

(a) 

(a) 

1“ 

0.005 (0.002) 

0.425 (0.008) 

1.7 (0.3) 

0.43 (0.08) 

3 (3) 

2° 

0.000 (0.001) 

0.392 (0.005) 

1.4 (0.2) 

0.33 (0.03) 

1.1 (0.9) 

3“ 

-0.003 (0.003) 

0.44 (0.01) 

0.7 (0.3) 

0.2 (0.1) 

1(3) 

Mean 0 

0.0013 (0.0009) 

0.407 (0.004) 

1.4 (0.1) 

0.36 (0.03) 

1.0 (0.8) 

4 ft 

0.022 (0.002) 

0.52 (0.01) 

1.8 (0.3) 

0.52 (0.02) c 

2.0 (0.5) c 

8 b 

0.023 (0.004) 

0.41 (0.02) 

1.9 (0.6) 

0.39 (0.03) c 

3 (1)° 

Mean 6 

0.021 (0.002) 

0.499 (0.009) 

1.9 (0.3) 

0.49 (0.02) 

2.5 (0.5) 


a Stowed-canard configuration. 
b Deployed-canard configuration. 
c Frequency domain estimate adjusted for bias. 
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Figure 18. Comparison of Z-force coefficient model parameters. Error bars represent la 
standard deviation. 
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4.5 Results Summary 

The stepwise regression method described in Section 3.2.3 reveals that linear models 
can adequately describe the aerodynamic pitch and yaw moment coefficients (Eqs. 19 and 
20, respectively), and body Y-force and Z-force coefficients (Eqs. 21 and 22, respectively). 
Models for the body X- force coefficient, and roll moment coefficient were not obtained. 


C m — C mQ + Cm a/ (a — 7 r) 

(19) 

C n = C no + Cnp/3 

(20) 

6V = 6V 0 +6V,/3 + 6V r ^ 

(21) 

Cz = C Zo + C Za , (a — 7r) + C Zq 

(22) 


Table 11. Parameter estimates obtained for the models for C m , C n , CV, and Cz for both 
stowed-canard and deployed-canard configurations. 


Parameter 

Time Domain 

Frequency Domain 

Stowed 

x(er) 

Deployed 

x(cr) 

Stowed 

x(ct) 

Deployed 

x(ct) 

Cm 0 

-0.006 (0.001) 

-0.069 (0.002) 

- 

- 

Cm a , 

-0.513 (0.007) 

-0.631 (0.008) 

-0.505 (0.006) 

-0.65 (0.01) 

Cn o 

0.009 (0.001) 

-0.006 (0.001) 

- 

- 

Cn, 

-0.561 (0.007) 

-0.585 (0.008) 

-0.532 (0.005) 

-0.59 (0.01) 

Cy 0 

0.001 (0.001) 

-0.002 (0.001) 

- 

- 

Cy, 

-0.422 (0.006) 

-0.487 (0.007) 

-0.46 (0.03) 

-0.48 (0.01) 

c Yr 

-2.2 (0.2) 

-1.5 (0.2) 

-3.5 (0.9) 

-1.5 (0.9) 

Cz 0 

0.001 (0.001) 

0.021 (0.002) 

- 

- 

Cz a , 

0.407 (0.004) 

0.499 (0.009) 

0.36 (0.03) 

0.49 (0.02) 

Cz q 

1.4 (0.1) 

1.9 (0.3) 

1.0 (0.8) 

2.5 (0.5) 


Parameter estimates for Eqs. 19-22 are provided in Table 11. A graphical comparison of 
these parameters can be found in Figure 19. The estimates and uncertainties are calculated 
using a covariance weighted average (see Section 3.2.4) of parameter estimates obtained 
for each test run. The estimates shown in Table 11 are displayed based on the significant 
digits of the uncertainties. Frequency domain estimates for the constant terms are not 
shown in Table 11. Large constant terms can bias frequency domain estimates toward low 
frequency content. This is prevented by removing these biases before transforming data to 
the frequency domain. The constant terms are not compared in Figure 19. 

Figure 19 compares the estimates for parameters determined in both the time and fre- 
quency domains. The close agreement between time and frequency domain estimates indi- 
cate that the model structures described by Eqs. 19-22 are accurate. Figure 19 also displays 
the effect that the canards have on the model parameters. Figure 19 and Table 11 show that 
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Figure 19. Graphical comparison of model parameter estimates in the time and frequency 
domains. Error bars represent la standard deviation. 
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the canards change the magnitude of all model parameters except for the constant terms, 
Cn,Q and Cy 0 • 

Table 12 provides each model’s Mean Squared Fit Error (MSFE) and R 2 for each test run. 
Table 12 shows that the identified models describe the pitch and yaw moment coefficients 
very well ( R 2 > 90%) with low fit error. The body Y-force and Z-force coefficients are not 
as well described. However, the MSFE for Cy and Cz are only slightly greater than the 
MSFE for C m and C n . 


Table 12. Mean Squared Fit Error (MSFE) and coefficient of determination ( R 2 ) describing 
the model fits for C m , C n , Cy, and Cz- 


Test 

Crr, 

l 

c n 


Cy 

Cz 

MSFE 

R z 

MSFE 

R z 

MSFE 

R z 

MSFE 

R z 

1“ 

0.003 

91.1 

0.001 

97.9 

0.002 

95.2 

0.001 

95.4 

2“ 

0.001 

97.9 

0.001 

97.6 

0.001 

95.9 

0.001 

96.7 

3 a 

0.002 

96.0 

0.001 

97.4 

0.004 

82.0 

0.003 

91.0 

3b b 

0.002 

94.2 

0.001 

96.7 

0.003 

84.9 

0.003 

85.7 

4 C 

0.002 

95.2 

0.001 

96.0 

0.0004 

98.0 

0.0005 

98.1 

8 C 

0.001 

95.7 

0.001 

93.4 

0.002 

80.7 

0.002 

81.2 

7 d 

0.001 

95.5 

0.001 

95.8 

0.002 

91.6 

0.001 

94.0 


“Stowed-Canard Configuration Analysis Data 
^Stowed-Canard Prediction Test 
c Deployed-Canard Configuration Analysis Data 
^Deployed-Canard Prediction Test 


Figures 20-23 display time histories and spectral densities of actual data, model predic- 
tions, and residuals for C m , C n , Cy, and Cz for the stowed-canard and deployed-canard 
configurations. The data used in Figures 20-23 was collected from test runs that were not 
used to determine parameter estimates (Tests 3b and 7 in Table 12). 

Figures 20 and 22 show how well the C m and C„ models fit the data as indicated in 
Table 12 for the stowed and deployed canard configurations. Notice that the residuals have 
higher frequency and much smaller values in the time domain plots. The spectral densities 
show that the models account for the main frequency content. The spectral density of the 
residuals is small and nearly flat indicating that the remaining residuals may be just noise. 

Figures 21 and 23 display time histories and spectral densities of the Cy and Cz for 
the stowed and deployed canard predication cases. As indicated in Table 12, the model fits 
for Cy and Cz are not as good as the fits for for the pitch and yaw moment coefficients 
discussed earlier. The larger residuals seen in the time history plots may indicate missing 
signal content in the data. The spectral densities plot show that models still account for 
the main frequency content of the actual data. 
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Figure 20. Moment coefficient predications for the stowed-canards configuration. 
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Figure 21. Force coefficient predictions for the stowed-canard configuration. 
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Figure 22. Moment coefficient predications for the deployed-canards configuration. 
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Figure 23. Force coefficient predictions for the deployed-canard configuration. 
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5 Summary 

A free-flight experiment concerning the Launch Abort Vehicle was conducted in the 
NASA Langley Research Center 20-ft Vertical Spin Tunnel. The goal of the experiment is 
to understand the aerodynamic characteristics of the Launch Abort Vehicle in the subsonic, 
incompressible flow regime. This regime represents a portion of the flight trajectory that the 
Launch Abort Vehicle experiences during the pad abort scenario. Two model configurations 
were flown during the experiment to determine how canards would impact the aerodynamic 
characteristics of the LAV. 

The primary objective of this experiment is the modeling of the aerodynamic pitch and 
yaw moments to determine static and dynamic stability of the Launch Abort Vehicle. Models 
were obtained for the pitch and yaw moments which indicate that the LAV is statically stable 
longitudinally and directionally in the heat shield forward flight state. Data analysis was 
unable to determine the longitudinal or directional dynamic stability parameters of the 
LAV. Models obtained for the canards show a change in trim angle to 174 degrees from 180 
degrees when no canards are deployed. 

Secondary objectives for the experiment included modeling aerodynamic forces and roll 
moment acting on the Launch Abort Vehicle. Models were obtained for the force coefficients 
acting in the body y- and z-directions. Models were not identified for the body x-direction 
force coefficient and the rolling moment coefficient. 
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